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p^ : Abstract 

I We propose a method for the stabihsation of quantum computations (including quan- 

^ ' tum state storage). The method is based on the operation of projection into SyA4, 
the symmetric subspace of the full state space of R redundant copies of the computer. 

^ ! We describe an efficient algorithm and quantum network effecting iSJ^A^-projection 

00 '. cind discuss the stabilising effect of the proposed method in the context of unitary er- 

g . j,Qj,g generated by hardware imprecision, and nonunitary errors arising from external 

' environmental interaction. Finally, limitations of the method are discussed. 

o : 
\o . 

: §1 Introduction 

p { 

^ I Any realistic model of computation must conform to certain requirements imposed 

Ch ' not by the mathematical properties of the model but by the laws of physics. Compu- 

^ ' tations which require an exponentially increasing precision or exponential amount of 

O^. time, space, energy or any other physical resource are normally regarded as unrealis- 

*> ! tic and intractable. 



^ I Any actual computational process is subject to unavoidable hardware imprecision 

I and spurious interaction with the environment, whose nature is dictated by the laws 

of physics. These effects introduce errors and destabilise the progress of the desired 
computation. It is therefore essential to have some method of stabilising the compu- 
tation against these effects. 

For classical computation there is a simple and highly effective method of stabili- 
sation. Each computational variable is represented redundantly using many more 
physical degrees of freedom than are logically required, and a majority vote (or av- 
erage) of all the copies is taken followed by resetting all the copies to the majority 
answer. This process is applied periodically during the course of the computation. If 
we use R copies and the probability of producing the correct answer is ^ + 1] then it 
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can be shown page 258) that the probabihty E that the majority vote is wrong, 
is less than exp(— r^^i^/G). This is an extremely resource-efficient stabilisation in that 
the probability of error decreases exponentially with the degree of redundancy R. 
Indeed, suppose that a polynomial-time algorithm runs for M steps, each of which 
is correct with probability \ + f] and majority voting is used after each step. The 
probability that the final answer will be correct is greater than (1 — E)'^^ . Thus 
any desired success probability 1 — 5 may be achieved using a degree of redundancy 
R = 0{log{M/6)), which is only logarithmic in the input size. 

The majority vote method above cannot be applied in the case of quantum computa- 
tion because quantum algorithms depend essentially on the maintenance of coherent 
superpositions of different computational states at each step. The laws of quantum 



mechanics forbid the identification of an unknown quantum state [IT^ [jT8[ and forbid 



even the cloning of an unknown state |jT9l. Thus the majority voting method is inap- 
plicable as we can neither determine the majority state nor reset the remaining copies 
to that state. In this paper we propose an alternative entirely quantum mechanical 
method of stabilisation utilising redundancy but which has no classical analogue. We 
discuss its applicability and limitations. The method was ffist proposed by Deutsch 
and a brief outline of its underlying principles was given in [Q]. 

An alternative approach to the stabilisation of classical computation involves the use 
of error correcting codes ^0\. A quantum mechanical generalisation of this approach 



was recently introduced by Shor and subsequently developed in [|T3l [|TOl |T2[ ||TT] . 
These methods are unrelated to those proposed in this paper and provide an inter- 
esting supplementary method of stabilising quantum computation. 

The process of simply repeating a whole computation a sufficient number of times 
may serve to stabilise it in certain circumstances. Suppose that we have a quantum 
algorithm which succeeds with probability 1 — e (where e may increase with input size) 
and suppose that we know when the computation has been successful. For example 
the computation may produce a candidate factor of an input integer which can then 
be efficiently checked by trial division. For any input size L the success probability 
can be amplified to any prescribed level 1 — 5 by repeating the algorithm a sufficiently 
large number, R, of times since the probability of at least one success in R repetitions 
is 1 — — 1 as -R ^ oo. Suppose now that the success probability 1 — e decreases 
with input size L as l/poly(L). Then we can maintain any prescribed success prob- 
ability by allowing R to increase as a suitable polynomial function of L. Thus if the 
original algorithm was efficient (i.e. polynomial time) then its -R-fold repetition is 
still efficient i.e. the algorithm has been stabilised in an efficient manner. (Note that 
Shor's quantum factoring algorithm is of this type with success probability 
decreasing as 1/L with input size.) However if the success probability falls exponen- 
tially with input size L then we must use R ~ exp(L) to maintain any constant level 
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of success probability. This implies an exponential increase of physical resources for 
stabilisation and hence this method is inefficient in this case. Unfortunately just such 
an exponential decay of success probability appears to be a generic feature of any 
physical implementation of computation, as described below. 

In quantum theory, the issue of preventing information from leaking into the environ- 
ment from a system ( "the computer" ) is generally known as the decoherence problem 
0, H, ^ . According to the analysis of , decoherence generally causes an exponential 
decrease in success probability with input size L. Decoherence is a universal phe- 
nomenon and is expected to affect - to some extent - any physical implementation of 
quantum computation whatever. Thus without some efficient form of stabilisation, 
quantum algorithms which are polynomially efficient in the error-free case (like Shor's 
factoring algorithm) cannot be considered polynomially efficient in practice. 

Consider any efficient computation which gives the correct result at each step with 
probability 1 — e where e is constant. This would typically be the case if each step 
of a computation consisted of the application of an elementary gate operation hav- 
ing a standard tolerance of error. Then after N steps the probability of success is 
(at least) (1 — e)^ ~ exp(— eA^) which again decreases exponentially with A^. Sup- 
pose that we have a stabilisation scheme utilising redundancy R which reduces the 
error in each step only by a factor 1/R i.e. e ^ e/R (rather than the exponential 
decrease given by classical majority voting). After A^ steps the probability of suc- 
cess is now exp{—eN/R). This can be kept at any prescribed level 1 — 5 by taking 
R = eN/ (— log(l — 6)) which is polynomial in A^ and hence in the input size L. Thus 
an exponentially growing error (such as results from decoherence) in a polynomial- 
time computation can be efficiently stabilised by a method which reduces the error 
per step only as with the degree of redundancy. Our proposed method below 
will have this property. 

It will be useful in the following to keep in mind the simplest possible example of 
the stabilisation problem where the computer consists of one qubit [i.e. one two- 
level system) and is performing no computation. In fact this simple model captures 
the essential features of the stabilisation problem for general quantum computations. 
The problem of stabilisation concerns the time evolution of an "accuracy" observ- 
able which has only two eigenvalues. As we shall see our analysis of error correction 
depends only on such simple observables and is independent of the substance of the 
computation. Thus we are addressing the problem of stabilising the storage of an 
(unknown) quantum state of one qubit against environmental interaction and (suit- 
ably random) imprecision in the construction of the storage device. 

§2 The Symmetric Subspace 
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Our proposed stabilisation method will exploit redundancy but in contrast to the 
classical majority voting method, it will be based on essentially quantum mechanical 
properties through use of the symmetric subspace of the full state space of R copies 
of a physical system. Consider R copies of a quantum system each with state space 
Ti. Denote the full state space Ti ®T-L ® . . .Tihy H^. 

Remark 1. Here we require that the R copies be distinguishable e.g. being located in 
separate regions of space so that the position coordinate provides an extra "external" 
degree of freedom for distinguishability. The state space 7i can be thought of as 
representing the "internal" degrees of freedom of each system. In our application 
these are the computational degrees of freedom of each replica of the computer. In 
our notation we suppress explicit mention of the distinguishing degree of freedom 
which is implicitly given by the written order of component states in a tensor product 
state (c.f. remarA; ^ below). □ 

The symmetric subspace SyM. of Ti^ may be characterised by either of the two 
following equivalent definitions. 

Definition 1. SyM. is the smallest subspace of Ti^ containing all states of the form 
\ip)... \^) for all |V) eU. n 

Definition 2. SyM. is the subspace of all states in Ti^ which are symmetric {i.e. 
unchanged) under the interchange of states for any pair of positions in the tensor 
product. (Here we interchange only the internal degrees of freedom leaving the external 
degrees fixed.) □ 

Remark 2. To clarify the notion of symmetrisation in definition 2 note that, for 
example, |0) + \ip) |0) G is in Sy^A. If we were to show the external degrees of 
freedom then this state would be written |0; xi) X2) + IV'; ^i) 10; X2)- Consequently 
the notion of symmetrisation in definition 2 is different from bosonic symmetrisation 
which requires symmetrisation of all degrees of freedom. For a pair of bosons the 
previous state would be |0; xi) X2) + I'?/'; X2) |0; xi). □ 

Definition 1 has the following interpretation. Suppose that we have R copies of a 
quantum computer. If there were no errors then at each time the joint state would 
have the form 

\^)\^)...\^)en'' (1) 

In the presence of errors the states will evolve differently resulting in a joint state 
of the form |?/'2) • • • iV'i?) or more generally a mixture of superpositions of such 
states. In quantum mechanics any test ("yes/no" question) that we can apply to a 
physical system must correspond to a subspace of the total state space. States of the 
form (JID for all \tp) eH do not, by themselves, form a subspace of H^. According to 
definition 1, Syj^ is the smallest subspace containing all possible error-free states. 
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It thus corresponds to the "most probing" test we can legitimately apply, which will 
be passed by all error-free states. Recall that we cannot generally identify the actual 
quantum state during the course of the computation or indeed gain any information 



about it without causing some irreparable disturbance |]I8[. The characterisation 
given in definition 2 is especially useful in treating mathematical properties of SyA4 
as below. 

The equivalence of the two definitions may be proved by viewing the i*^ component 
space in the tensor product TC^ as the space of complex polynomials of degree < n — 1 
in the variable (where n is the dimension of Ti). Then Definition 1 defines the 
subspace of all polynomials p{xi, . . . , xr) of degree < n — 1 m. each variable, which 
arise as sums of products of functions of the form 

R 

fa{xi,X2,. . . ,Xr) = \{{xi- a) (2) 

i=l 

for any a (where we have used the fundamental theorem of algebra). On the other 
hand, definition 2 defines the space of all symmetric polynomials (of degree < n — 1 
in each variable). The equivalence of these subspaces then follows easily from basic 
properties of the standard elementary symmetric functions Jill which are defined as 
the coefficients of the powers of a in the expansion of (^. 

The equivalence of the two definitions may also be understood via the following illus- 
trative example which gives further insight into the structure of SyM.. 

Example 1. Suppose that Ti is two-dimensional [i.e. a qubit) with computational 
basis states |0) and |1). Consider triple redundancy R = 3 and the symmetric sub- 
space Syj^ C H^. Let us tentatively denote the symmetric subspaces of definitions 1 
and 2 by Syj^defi and SyM.def2 respectively. We wish to show that these coincide. 
Note first that Sy^Adefl is the span of all states of the form \ip) lip) \ip), which are 
clearly symmetrical in the sense of definition 2. Hence Syj^defi ^ Syj^def2- For 
the reverse inclusion consider a general state in Ti^: 

\a) = ao|0)|0)|0) 

+ai |1) |0) |0) + a2 |0) |1) |0) + as |0) |0) |1) (3) 
-Fa4|l)|l)|0) + a5|l)|0)|l) + a6|0)|l)|l) 
+«7|1)|1)|1) 

Interchange of states for any given pair of positions (in the sense of definition 2) 
preserves the number of |0)'s and |l)'s in each term so that \a) will be in SyAidef2 
if and only if ai = 02 = and = = a^. Indeed we see that SyA4def2 is four 
dimensional with orthonormal basis states (labelled by the number of |l)'s): 

|eo) = |0)|0)|0) 
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|ei) = (|l)|0)|0) + |0)|l)|0) + |0)|0)|l))/v^ 
|e2) = (|1)|1)|0) + |1)|0)|1) + |0)|1)|1))/V3 
lea) = |1)|1)|1) 



(4) 



(The four normalising factors 1, -\/3, and 1 are square roots of the binomial coef- 
' Ci,^ 6*2,^ C3.) Now for any of the form a |0) + |1) we get directly 



ficients 
that 



0; 



li^i) l"^!) iV^i) = «^ |eo) + a^Vs |ei) + aVs {62) + {e^) 



Repeating this for four different values of the parameter a we get: 



/ a3 



a 



a 
b 
c 
d 



1\ 

1 

1 
1 / 



f |eo) \ 

73162) 

V 1^3) J 



( 1^1) \ 
1^2) 1^2) 1^2) 

iV^s) 1^3) 1^3) 

V iV'4) 1^4) 1^4) ; 



Choosing a, 6, c and d so that the coefficient matrix is invertible, we see that the basis 
states are all in Syj^defi so that SyM.def2 ^ SyM.defi- Hence these subspaces 
coincide. □ 



From the above considerations (c.f. especially (D ) we readily see that the dimension 
of SyM. for R qubits is i? + 1 so that SyM. is an exponentially small subspace 
of Ti^ (of dimension 2'^). This is also true in the general case. Suppose that Ti has 
dimension d with orthonormal basis |0) , |1) . . . |(i — 1). Then SyM. has an orthogonal 
basis labelled by all possible ways of making R choices from the d basis states with 
repetitions possible and the ordering of choices being irrelevant (c.f. (|^) and @)). 
The solution of this combinatorial problem gives 

Dimension of SyM = ^^'^'^Ca-i = tt^-^R"^'^ + 0(R'^-^) (5) 

{d-iy. 

which is a polynomial in R (for fixed d). Hence SyAi is again exponentially small 
inside Ti^ of dimension d^. 



§3 Projection into SyA4 

Our proposed method of stabilisation consists of frequently repeated projection of 
the joint state of R computers into the symmetric subspace SyAi. According to 
the interpretation of SyA4 above, the error free component of any state always lies 
in Syj^ so that upon successful projection this component will be unchanged and 
part of the error will have been removed. Note however that the projected state is 
generally not error-free since, for example, SyAi contains many states which are not 
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of the simple product form \ip) \ip) . . . I^^). Nevertheless the error probability will be 
suppressed by a factor of l/i? as discussed in subsequent sections. Thus the method 
is not one of error correction but rather of stabilisation. By choosing R sufficiently 
large and the rate of symmetric projection sufficiently high, the residual error at the 
end of a computation can, in principle, be controlled to lie within any desired small 
tolerance. 

The operation of projection into SyM. is a computation in itself. For our stabilisation 
method to be efficient it is essential that this operation be executable efficiently i.e. in 
a number of steps that increases at most polynomially with L and R where L = logg d 
is the number of qubits required to hold the state of each computer entering into the 
symmetrisation and R is the degree of redundancy. (Note also that clearly R can 
be at most a polynomial function of L in any efficient scheme.) Only then will a 
nominally efficient computation remain efficient after stabilisation. 

We next describe an algorithm for projecting into SyAi and show that it is efficient 
in the above sense. Consider first a product state |^) = |ai) \a2) ■ ■ ■ Iqr) G 7i^. To 
project 1^) into SyAi we carry out the following steps: 

Step 1: Introduce an ancilla in a standard state |0) with a state space A of at least 
7?! dimensions. 

Step 2: Make an equal amplitude superposition of the ancilla 

1 m-i 

^ HO) - ^ E N) 

Step 3: Carry out the following computation: if the ancilla state is \i) then perform 
the i*'* permutation Uj of the component states of |ai) [02) . . . {an) 



\ai) \a2) . . . \aR) \i) - 
This results in the entangled state 

E |«'^i(2)) • • • \i) en^^A 

i 

Step 4: Apply the reverse computation of step 2 to the ancilla. The resulting 
state may be written 

\r)^Y:\^i)\i)en^0A 

i 

Since U transforms |0) to each \i) with equal amplitude it follows that U^^ 
transforms each \i) back to |0) with equal amplitude. Hence the coefficient of 
ancilla state |0) in |T) is the required symmetrised state i.e. an equal amplitude 
superposition of all permutations of the R factor states of |ai) |a2) . . . 
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Step 5: Measure the ancilla in its natural basis. If the outcome is "0" then 

has been successfully projected into SyA4. If the outcome is not "0" then the 
symmetrisation has failed. (The issue of the probability of successful symmetri- 
sation is discussed in a later section.) 

Finally note that by linearity of the process, it will symmetrise a general state in Ti^ 
(not just the product states considered above). If the input state is already symmetric 
then we get it back unchanged with certainty at the end. 

Consider now the computational effort involved in the above steps. Let d = dim 7i 
and write L = log2(i. Step 1 requires no computational effort. The ancilla requires 
log2(-R!) = 0{RlogR) qubits. Step 2 may be achieved by applying the discrete 
Fourier transform HI to the ancilla. This requires 0((-Rlog/2)^) steps. For step 3 
we note that a general permutation can be effected with 0{R\ogR) swaps. Swapping 
states of L qubits requires 0{L) operations so overall step 3 requires 0{LR\ogR) 
steps. Restoring the ancilla in step 4 requires the same number of operations as step 
2. In step 5 we examine separately each of the 0{R\ogR) qubits occupied by the 
ancilla, requiring 0{R\ogR) steps. Overall we require 0(L_Rlog-R+ (-Rlog/2)^) steps 
which is less than 0{LR^ + i?^). Hence the process is efficient. 



§4 A Quantum Network for SyA4 Projection 

We now describe how the operation of SyAi projection can be implemented by a 
network of simple quantum gates. Consider first the following inductive definition 
of the general permutation of A; + 1 elements ai, . . . , a^, afc+i[|^. Starting from the 



general permutation ao-(i), . . . , a^(^k) of the k elements Oi, . . . , we adjoin a^+i giving 
ao-(i) , . . . , a^{k) , flfc+i and then perform separately the k + 1 operations: identity, swap 
ao-(i) with ttk+i, swap a^(^2) with a^+i, . . . swap ao-(k) with ak+i- This generates all 
possible permutations of + 1 elements. In terms of state symmetrisation, if we have 
already symmetrised \ipi) (S) . . . ® li'k) (^-e. we have an equal superposition of all per- 
mutations of the states) then we can symmetrise k + 1 states \ipi) Cg) . . . CS) \ipk) ® |^fc+i) 
by applying only the operation of state swapping (in suitable superposition). Thus 
the operation of symmetrisation of R states can be built up from by first sym- 
metrising \ipi) and \4'2), then successively including lip^) up to lipji) always using only 
state swappings in suitably controlled superpositions. 

The basic ingredient in this process is the "controlled swap gate" or Fredkin gate, 
acting on three input qubits. If the first ("control") qubit is |0) (respectively |1)) 
then the other two ("target") qubits are unaffected (respectively swapped). We de- 
scribe this diagramatically in Fig. 0. The operation of state swapping itself {i.e. 
I'ipi) ® IV'2) ^ IV'2) ® 1"^!)) can be implemented using three controlled-NOT gates as 



described in [15 
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Figure 1: Schematic representation of a Fredkin gate. A Fredkin gate exchanges the state 
of the second and third qubit if and only if the first qubit is in state \a) = 



To symmetrise k + 1 qubits given that the first k are already symmetrised we 
introduce k control qubits initially in state |0) |0) . . . |0) and apply a suitable unitary 
transformation, denoted Uk to generate the superposition 



1 



(|00 . . . 0) + |10 . . . 0) + |01 . . . 0) + • • • + |00 . . . 1)) 



(6) 



The unitary transformation Uk can be readily obtained by a quantum network con- 
sisting of a one bit gate performing the transformation 



1 ( 1 -Vk\ 
Vk + T \Vk 1 ) 

on the first qubit and a sequence of — 1 two bit gates Tjj^i for j = 1, . . 
on the and j + 1*'* qubits. In the basis {|0) , |1)}, Tjj+i is given by: 



(7) 

k — 1 acting 



T 1 



( Vk-j + 1 



V 





1 










■Vk^ 1 

^/k-j + 1 J 



(8) 



Having thus initialised the k control qubits, we then apply k Fredkin gates: the j*^ 
Fredkin gate (for j = 1, . . .k) uses the j*^ control qubit to control the swapping of 
the j*'* and k + 1*'* target qubits. This leads to an entangled state of the k control 
qubits and the A; + 1 target qubits but after applying U^^ to the control qubits, the 
coefficient of |0) |0) . . . |0) will be the required symmetrisation of the k + 1 qubits (c.f. 
step 4 of §3). Finally a measurement of the control qubits will effect the projection 
into SyM (c.f. step 5 of §3). 
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Thus to symmetrise R qubits we cascade the above construction with k = 1,2... 
up to k = R — 1 requiring a total number l + 2 + -- - + (i? — 1) = R{R — l)/2 
of control qubits. The size of the overall network is clearly quadratic in R. For ex- 
ample, for the symmetrisation of i? = 4 qubits we obtain the network shown in Fig. 0. 



|0> 

lo) 
lo) 

lo) 
lo) 



U3 



lo) -Ui 



k) 
k) 
k) 
k) 



—7^ 

^ ^ 



^ ^ N/ 



U3 



U2 



/ 



Ui 



Figure 2: Quantum network for symmetrising ii = 4 qubits. Six auxiliary qubits initially in 
state |0) are needed. The auxiliary qubits are put into an entangled state and used to control 
the state swapping of the four computer qubits. The operations are then undone and the 
auxiliary qubits measured. If every auxiliary qubit is found in state |0) the symmetrisation 
has been successful. 

§5 Stabilisation Against Unitary Errors 

So far we have given an efficient algorithm for projection into the symmetric sub- 
space and provided an intuitive reason why it would be expected to reduce the error 
while preserving the correct computation. We now turn to a quantitative study of 
the effect of iS^A^-projection as a basis for stabilisation in the presence of various 
modes of error production. It is convenient to separate the discussion into two parts 
considering the case where the joint state of the computers remains in a pure state, 
in this section, and the case of decoherence due to external environmental interaction 
in the following section. 
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Consider the simple model of R qubits initially in state (the "correct" state) |0) |0) . . . |0) 
with computation being the identity i.e. we are considering the state storage prob- 
lem with _R-fold redundancy. Suppose that the R storage devices are subject to 
independent hardware errors which cause the j^^ state to drift as e*^^* |0). Here the 
Hamiltonians Hj are random and independent. Since the devices were intended for 
state storage we assume that the rate of drift is suitably bounded. This is expressed 
by requiring that all eigenvalues of the -ffj's are suitably small: 

jeigenvalues oi Hj\ < e j = 1, . . . , R (9) 

for some (small) constant e. The stabilisation process consists of projecting the joint 
state of the R copies into the symmetric subspace at short time intervals 6t. For 
simplicity we will assume that the projection can be performed essentially instanta- 
neously. Furthermore we assume that (unlike the computation being stabilised) the 
projection process itself is error free. These and other assumptions will be discussed 
in §7 below. Under these assumptions we can readily compare the growth of errors 
with and without the stabilisation process. 



In the basis {|0) , |1)} write 



so that 

|cj| < |Ai| + IA2I < 2e (11) 

(where Ai, A2 are the eigenvalues of Hj.) We will assume that 5t is small and retain 
only the lowest order terms in 6t. After time 6t the state will be 

R 

\^(St)) = (g) {(1 + ^akSt) |0) + zck6t |1)} (12) 

k=l 

Thus without symmetrisation the probability that the i*^ qubit shows an error is 

\c,\^5t^ ^ Ae^dt^ (13) 
If we expand out the product in (p!2|) we obtain 2^ terms corresponding to the expo- 



nentially large dimension of the full space of R qubits. However the amplitudes of 
terms involving k errors {i.e. products of |0)'s and |l)'s involving k |l)'s) will have size 
0{6t'') and only the R terms involving one error will have size 0{6t). Thus the erro- 
neous state ( p!2D does not occupy these exponentially many dimensions of equally. 
We noted previously that Syj^ is exponentially small inside TC^ but the preceeding 
observation indicates that iS^A^-projection will not generally remove exponentially 
much of the error since only R of these exponentially many dimensions are entered 
(to lowest order) by the erroneous evolution. We now calculate the stabilising effect 
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of the iS3^7Vl -projection. 

Consider the basis of SyAi given by the R + 1 orthonormal states (c.f. @ for the 
case R = 3): 

|e.) = ^^ E 1^) k = 0,...,R (14) 

^ ^'^ all "A;-error" a's 

Here the sum is over all ^C^ possible strings of O's and I's of length R containing 
exactly I's and R — k O's. Under -projection the lowest order error terms 

(single-error terms) of ([l^) will project only onto |ei). For the term with an error in 
the k^^ place we get 

• ;r+ln\ ln\h\ln\ ln\ 53^A^ proj iCkSt 
zcfcdt |0) . . . |0) |1) |0) . . . |0) ^ ^==|ei) 



Thus the normalised projected state has the form 

l+i6tJ2 % [ |0) • • • |0) + ai |ei) + 0{6t^) 



where 



Rj-_ 



To estimate the size of ai, using (|T^) we write Cj ~ 2ee*^^ where 6j are random phases. 
The expectation value of ai is then clearly zero but from 

R 

Expectation value of | E I = (16) 

i=i 

we get Q 

Expectation value of |ai|^ = Ae^St^ (17) 

Thus, somewhat surprisingly, this probability of a single symmetrised error does not 
decrease with R. However it is associated with R copies and to see its residual effect 
on any one copy we use the following fact: 



^ Remark : If instead of qubit systems we consider computers with dimensions large compared 
to the degree of redundancy R, then we would expect the individual random errors to be mutually 
orthogonal, so when the state is symmetrised their sum does not exhibit the cancelling effects 
which are present for qubits and lead to (|l^) and (|l^). However in that case, and subsequent 
conclusions still hold because ( p^ ) may be replaced by Pythagoras' theorem i.e. that the sum of R 
orthonormal vectors has length ^/R. 
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Proposition 1 Consider the state 



R 

\'E) = Y.au\eu)eSyM<zn'' 

k=0 

where \ek) are as given in (|14|). If one qubit is measured in the basis {|0) , |1)} the 
probabihty of obtaining |1) is ^J2k=o l^^feP^- '-' 

Proof: Since the state is symmetric the probabihty of obtaining the result |1) for the 
i^^ qubit is the same as this probabihty for the first qubit. Now \ek) in (|I^) consists 
of ^Ck orthogonal terms of which ^~^Ck~i have |1) in the first place. Hence the term 
«fc |efc) in |S) contributes probability 



of obtaining outcome 

Applying this result to (^) and using (|T7p we see that after successful symmetrisation 
the probability of error (to lowest order in 6t) is 4e'^5t^/R i.e. the error is suppressed 
by a factor oi 1/R compared to the case (|l^) of no symmetrisation and in each step 
the amplitude of correct computation is correspondingly enhanced. 

The above result is conditional on the success of the symmetrisation i.e. that the 
state projects to SyAi rather than SyAi^. If the projections are done frequently 
enough then the cumulative probability that they all succeed can be made as close 
as desired to unity. This is a consequence of the so-called "quantum watch-dog 
effect". Consider a normalised joint state |S) of R copies initially in SyA4. Its 
initial probability of successful projection is 1 which is a maximum. Thus as the state 
evolves by some unitary transformation into the ambient space Ti.^ the probability of 
successful projection will begin to change only to second order in time. If we project 
n times per unit time interval i.e. 6t = 1/n then the cumulative probability that all 
projections in one unit time interval succeed, is 

(1 - Mt^)- = (1 _ 4-)" ^ 1 as n ^ oo 

Here k is a constant depending on the rate of rotation of the state out of Syj^. For 
redundancy degree R and the model of random unitary errors considered above we 
find that k grows linearly with R (as can be seen by directly calculating the length of 
the 53^A4-projection of (O) to 0(5t^) terms). Thus to achieve a cumulative proba- 



bility of successful projection of 1 — C in a unit time interval we would require a rate 
of symmetric projection which increases linearly with — -R/ log(l — Q. 

The above conclusions — for a model of random independent unitary errors — will 
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also apply to computations which are not the identity. Formally we may view the 
computation in a moving basis relative to which the correct computation is the iden- 
tity and the previous arguments are unchanged i.e. none of the arguments depend on 
the actual identity of the computational basis states. 

§6 Stabilisation Against Environmental Interaction 

We now consider the problem of state storage with i?-fold redundancy, in the pres- 
ence of decoherence i.e. interaction with an external environment. In general each 
qubit will become entangled with an environment and the state of the qubit alone will 
no longer be describable by a pure state. It will be represented by a density matrix 
IE] resulting from forming a partial trace over the environment, of the joint (pure) 



state of the total qubit-environment system. 

Consider R copies of the qubit initially all prepared in pure state po = |0) (0|. We will 
assume that they interact with independent environments (this assumption is valid if 
the coherence length of the reservoir is less than the spatial separation between the 
copies so that after some short period of time 6t the state of the R copies will 
have undergone an evolution 

p(^)(0) = Po ® . . . ® Po ^ P^'^HSt) = Pi ® . . . ® PR (18) 

where Pi = Po + o"t for some Hermitian traceless (jj and the superscript R denotes the 
number of states involved. We will retain only terms of first order in the perturbations 
ai so that the overall state at time 6t is 

P^^^ = Po ® • • • ® Po +<^i ® Po ® • • • ® Po 

+P0 ® ® . . . ® Po (19) 

+Po ® Po ® ■ ■ ■ ® ctr 
+0{aiaj) 

and we wish to compute the projection of the state (|1^) into the symmetric subspace 
SyA4. Then we construct the state of the i^'^ qubit by partial trace over all qubits 
except the i^^ and finally compare the resulting state with po + cxj and see that its 
purity has been suitably enhanced, bringing it closer to po. 

The mathematical formalism for symmetrisation of mixed states has some curious 
features which we digress to clarify before treating (|r^) itself. Consider a state pi® P2 
of two qubits. The state \{pi ® P2 + P2 ® Pi) is not a symmetric state and in fact 
p ® p is not symmetric {i.e. it is not a density matrix supported on the subspace 
SyM) unless p is pure! To see this consider p written in its diagonalising basis of 
orthonormal eigenstates: 

P = Ai|Ai)(Ai| + A2|A2)(A2| (20) 
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Thus we can represent p as a mixture of its two eigenstates, and p p as a mixture 
of the four orthonormal states |Aj) ® \Xj) with a priori probabihties pij = XiXj. This 
latter mixture involves nonsymmetric states (like |Ai) ® IA2)) so p®p is not symmet- 
ric. 

One way of constructing the projection of p ® p into Syj^ is to project each state of 
the above mixture into SyA4. Let |pjj) denote the iSJ^A^-projection of |Aj) ® \Xj) 
and \fiij) denote the corresponding normalised state. The probability of successful 
projection is qij = (fJ^ijlfJ^ij)- Then the iS3^7Vl-projection of p ® p is the state corre- 
sponding to the mixture |pjj) with a priori probabilities Pijqij/ (J2Pij1ij) which are the 
conditional probabilities of occurrence of states |pjj) given that the iSJ^A^-projection 
was succcessful. 

More formally we may introduce the (Hermitian) permutation operators P12 = 
"identity" and P21 = "swap" acting on pure states of two qubits and define the 
symmetrisation operator: 

S=^{P,, + P,,) (21) 

The 53^A^-projection of a pure state |\E'i2) of two qubits is just 5* l^^u), which is then 
renormalised to unity. It follows that the induced map on mixed states of two qubits 
(including renormalisation) is: 

The state of either qubit separately is obtained by partial trace over the other qubit. 

As an example consider the symmetric projection of p®p followed by renormalisation 
and partial trace (over either qubit) to obtain the final state p of one qubit, given 
that the iS3^A^-projection was successful. A direct calculation based on (P^) yields: 



For any mixed state of a qubit the expression Tr S,'^ provides a measure of the purity 
of the state, ranging from ^ for the completely mixed state 1/2 (where I is the unit 
operator) to 1 for any pure state. From (53) we get 



Tr p > Tr p 

so that p is purer than p. This example illustrates a generic fact (c.f. below), that 
successful projection of a mixed state into SyM. tends to enhance the purity of the in- 
dividual systems. Indeed, consider further the state consisting of R independent 
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copies of p. The symmetrisation operator is 



5=^|:P. (24) 
where the sum ranges over all R\ permutations of the R indices. If we project 



into SyM. and renormalise (as in (]22|) ) and calculate the partial trace over all but 
one of the qubits, we obtain a reduced state pr which asymptotically tends to a pure 
state as R oo. This limiting pure state is the eigenstate of p belonging to its largest 
eigenvalue. 

Let us now return to the consideration of ([T9|) and its 53^A^-projection. The applica- 
tion of the symmetrisation operator (p^) to each of the R terms of po®- • .®(Tj®. . .®Pq 
of Eq. ([T9|) generates i?!^ terms of the form 

PapO ® . . . ® CXi ® . . . ® poPp, (25) 

where Pa and Pp are permutations operators on the state space Ti^ of R qubits as 
above. To calculate the reduced density operator of the first qubit we take the partial 
trace over the R — 1 remaining qubits. Note that the reduced states of all qubits 
individually are equal since the total overall state is symmetric. (To systematise 
the calculation of the partial traces we have found it very convenient to use the 
diagrammatic notation for tensor operations introduced by Penrose in [^.) For each 
(Tj we find that the i?!^ terms in (pSj) then reduce to the following cases: 



(i) {R — 1)!^ terms each equal to (JijW? corresponding to all permutations Pa and 
P(i which place ai in the first position in (pSj). In this case the partial trace 
contracts out all the po terms leaving a coefficient of XjR}? (as the trace of any 
power of po is 1). 

(ii) [R — 1)!^(-R— 1)-R terms of the forms poaiPQ, PoO'i, <JiPo, or poTr(crjpo), each one 
divided by R\^. These correspond to all pairs of permutations which result in 
(Ji contracted onto po in all possible ways in the partial traces. 

(iii) {R — 1)!^(-R — 1) terms which result in ai being contracted onto itself in the 
partial traces. These terms are all zero since Tr cr, = 0. 

Note that each term in (ii) has trace given by Tr aip/ R\^ and each term in (i) has zero 
trace. Thus the resulting density operator, before normalisation, has a trace given by 

l + (i?-l)Tr(poa), (26) 
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where we have introduced a = After normahsation - which is calculated 

by dividing the sum of all terms in (i) and (ii) for alH = 1, . . . i? by ( p6D - the resulting 
symmetrised density operator p can be written 

p = [1 - (i? - l)Tr(po^)]po + -^^ 

+ {R-l)[Apoapo + B{poa + apo) + CpoTT{dpo)]+0{aia,) (27) 

where A, B and C depend on R and A + 2B + C = 1. 

If a general mixed state ^ of a qubit is measured in the basis {|0) , |1)} then the 
probability that the outcome is is given by (0|^|0) = Trpo^- This provides the 
success probability in our present model. Thus the average success probability before 
symmetrisation of the perturbed qubits is: 

4ETrpo(po + ^^) = 1 + Trpo^ (28) 

(Note that consequently Tr poO" is necessarily negative). y4/i(er symmetrisation, using 
( p7D we see that 

Trpop = l + ^Trpo^ (29) 
ri 

Hence the probability of error has again been reduced by a factor of R - exactly as 
found in the previous section. 

We can calculate the average purity of the R copies before symmetrisation by calcu- 
lating the average trace of the squared states: 

^ X] Tr((po + (T.f) = 1 + 2Tr(poa). (30) 
After symmetrisation each qubit has purity 

Tr(p2) = l + 2iTr(poa). (31) 

Since Trp^ is closer to 1 than (^OD, the resulting symmetrised system p is left in a 
purer state. Indeed it follows from (|29| ) that p approaches the unperturbed state po 
as R tends to infinity. 

§7 Limitations 

Error correction is itself a quantum computation. The above analysis has ignored the 
inevitable build up of errors in the computer performing that computation. Indeed 
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for the symmetrisation of R qubits the projection algorithm requires an ancilla of at 
least R\ dimensions i.e. 0{R\ogR) qubits (in fact 0{R^) in our exphcit network). 
Thus the correcting apparatus is shghtly larger than the total system being corrected 
so the error correction ought to be subject to a similar level of error as is present in 
the original system. In a situation where the redundancy degree R is small compared 
to the number L of qubits per computer, the correcting apparatus (still of O^R?) 
qubits) will be small compared to the size RL of the system being corrected. How- 
ever as seen in §5, the stabilisation of a linear computation on input size L requires 
redundancy degree R ^ L so that the correcting apparatus and the computer are 
again of comparable size. This means that each error correcting step introduces er- 
rors of a similar, or even greater, probability than those that it is correcting. This 
does not however render it ineffective. Consider the following illustrative example. A 
certain clock is accurate to one second per day. Each day at noon it is reset using a 
standard time signal, the resetting operation being accurate only to one minute i.e. 
sixty times worse than the error being corrected. Nevertheless after ten years the 
corrected clock will still be in error by at most one minute. If left uncorrected the 
error could be almost an hour. 

The main factor limiting the efficiency of our proposed method will be the frequency 
with which the error correcting operations can be physically performed. As noted at 
the end of §5, to achieve a cumulative probability 1 — 5 of repeated successful projec- 
tion in a unit time interval, the rate of symmetric projection must increase linearly 
with the degree of redundancy R. Also as noted in §1, the stabilisation of a computer 
with input size L, running for L steps, requires R to increase linearly with L. Hence 
we need the overall rate of symmetric projection to increase linearly with L even for 
a linear time computation. Thus, beyond a certain input size, each symmetrisation 
will have to be performed in a a time shorter than that needed to perform the el- 
ementary quantum gate operations. Since increasing the rate of computation by a 
factor k presumably requires resources exponential in k, our method would necessarily 
require exponential resources for sufficiently large L. This property is shared by all 
quantum error correction schemes that have been proposed to date. Hence quantum 
algorithms (such as Shor's factoring algorithm) which are polynomially efficient in 
the absence of errors, would not be efficient if physically implemented. We wish to 
stress that the traditional notion of efficiency (based on the distinction between poly- 
nomial and exponential growth) is an asymptotic notion referring to computations on 
unboundedly large inputs. This may not be appropriate in assessing the feasibility 
of particular computations in practice. For example, if a quantum computer could 
factorise a 1000 digit integer in a reasonable time it may still exceed the abilities of 
any classical computer for the foreseeable future albeit that the factorisation of 2000 
digit integers might be infeasible on any computer. 
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Conclusion 

If the technology to implement the scheme we have described were available, it would 
provide a method of stabilising general coherent computations though not (exponen- 
tially) efficiently. This is because although only polynomially many steps are required 
in the stabilisation computation, these need to be performed in a fixed time, a char- 
acteristic time of error growth per bit. 
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